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1.  INTRODUCTION 


Quasi-static  deformations  in  response  to  atmospheric  pressure  changes  account 
for  much  of  the  ambient  earth  noise  at  periods  greater  than  20  seconds.  While 
local  variations  in  atmospheric  pressure  may  arise  from  a variety  of  sources, 
only  two  appear  to  produce  significant  earth  motion  on  a more  or  less  continu- 
ous basis  These  are  the  turbulent  air  flow  associated  with  the  surface  wind 
(Sorrells  et  al.,1971,  Ziolkowski,  1973)  and  naturally  occurring  infrasomc 
waves  (Sorrells  and  Douze,  1974).  The  essential  difference  between  these  two 
phenomena,  insofar  as  the  earth  is  concerned,  is  that  they  occupy  widely 
separated  regions  of  the  wave-number  spectrum  of  the  atmospheric  pressure  “eld. 
In  the  20-100  second  period  range,  the  convective  wavelengths  associated  with 
wind  generated  turbulence  rarely  exceed  a kilometer  and  may  be  as  s„w  1 1 as 
several  tens  of  meters.  On  the  other  hand,  infrasonic  pressure  fluctuations 
in  the  same  period  range  commonly  have  wavelengths  which  are  of  the  order  o 
several  tens  of  kilometers.  Since  the  quasi-static  pressure  response  of  the 
earth  tends  to  increase  monotonically  as  the  wavelength  of  the  pressure  oscilla- 
tion increases  (Sorrells,  1971,  Sorrells  et  al.,1971),  infrasonic  waves  will, 
in  general,  cause  much  larger  deformations  than  wind-generated  atmospheric 
changes  with  the  same  power  spectral  density.  At  the  surface,  aunng  windy ^ 
intervals,  this  differencf  in  responses  is  more  than  offset  by  the  differences 
in  the  local  pressure  amplitudes  of  'he  two  sources.  Because  of  the  differences 
in  wavelengths,  the  earth's  response  to  acoustic  waves  may  exceed  its  response 
to  wind  related  pressure  changes  by  about  a factor  of  5-10.  However,  the 
amplitudes  of  wind  related  atmospheric  pressure  fluctuations  can  easily  be 
1-2  orders  of  magnitude  greater  than  the  amplitudes  of  infrasonic  waves  with 
the  same  period.  Therefore,  during  intervals  of  moderately  high  wind  speeds 
(>4  meters/sec)  the  pressu’ 2 related  earth  noise  recorded  at  or  near  the  sur- 
face will  be  dominated  bv  quasi-static  deformations  in  response  to  wind- 
generated pressure  changes.  Both  theoretical  studies  (Sorrells,  1971, 
Ziolkowski,  1973)  and  experimental  studies  (Sorrells  et  al ^emonstarate 
that  wind  related  earth  noise  can  be  vertu3lly  eliminated  in  the  20-100  second 
period  range  by  installing  the  seismographs  at  depths  on  the  order  of  several 

hundred  meters. 


An  alternative  method  of  reducing  pressure  related  earth  no<se  was  also 
considered.  It  has  been  known  for  some  cime  that  during  windy  intervals 
there  is  a strong  correlation  between  the  outputs  of  colocated  microbarographs 
and  seismographs  (Capon,  1969,  Sorrells  et  al . , 19/1).  Therefore,  it  is 
possible  to  estimate  the  pressure  generated  earth  noise  by  applying  an  appro- 
priate filter  to  the  output  of  the  microbarograph.  Noise  reduction  is  obtained 
by  subtracting  this  estimate  from  the  recorded  seismogram. 
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Information  collected  from  three  experimental  long-period  installations  will 
be  used  in  this  report.  A brief  description  of  each  site  follows: 

a.  McKinney,  Texas;  consisted  of  meteorological  equipment  and  a surface 
seismological  observatory,  and  was  used  principally  to  study  the  effect  of 
wind-induced  pressure  variations  on  surface  seismographs.  The  meteorological 
instrumentation  consisted  of  13  microbarographs  and  three  anemometers; 

one  microbarograph  array  of  six  instruments  was  closely  grouped  around  the 
seismometer  vault  with  average  separations  of  around  1 km  and  included  a 
mic robarograph  colocated  with  the  seismograph  vault.  The  remaining  seven 
microbarographs,  with  separations  of  around  5 km  between  sensors  constituted 
a second  array,  also  available  to  study  long-wavelength  acoustic  waves.  Four 
microbarographs  were  installed  for  a limited  time  at  a distance  of  50  m from 
the  vault  to  study  wind-generated  noise.  The  basic  seismic  instrumentation 
consisted  of  a pair  of  horizontal  and  a pair  of  vertical  seismographs, 

Teledyne  Geotech  Model  7505A  and  Geotech  solid-state  amplifiers  (Model  28450). 
These  surface  seismometers  were  enclosed  within  insulated,  pressure-tight 
tank  vaults  which  were  mounted  in  the  floor  of  a concrete  bunker.  The  roof 
of  the  bunker  was  below  ground  level  with  access  through  a pressure-tight 
hatch.  Data  from  all  the  sensors  were  recorded  on  a digital  acquisition 
system  described  in  detail  by  Herrin  and  McDonald  (1971).  In  addition  the 
data  were  recorded  continuously  on  16  mm  film. 

b.  Grand  Saline,  Texas;  experiments  were  performed  at  a salt  mine  owned 
by  the  Morton  Salt  Company.  The  seismic  instrumentation  consisted  of  two 
3-component  sets  of  long-period  seismographs  installed  at  the  surface  and  in 
the  mine.  The  mine  seismographs  were  at  a depth  of  183  m below  the  surface 
and  approximately  100  m of  horizontal  distance  from  the  surface  installation. 

One  microbarograph  was  colocated  with  the  surface  system  and  one  microbarograph 
was  used  to  measure  pressure  variation  inside  the  mine.  The  system  is  described 
in  detail  by  Herrin  and  McDonald,  1971,  and  Sorrells  et  al . , 1971. 

Directly  below  the  site  the  velocity  structure  is  controlled  by  the  salt  dome. 
Approximately  90  m of  sandy  clay  with  interbedded  sand  lenses  occur  above  the 
salt.  Compressional  wave  velocities  of  1.4  and  4.7  km/sec  (Anderson  and 
Lieberman,  1946)  were  used  for  the  salt  dome  mode  with  a Poisson's  Ratio  of  0.33. 
The  velocity  structure  of  the  sediment  surrounding  the  salt  dome  were  taken  from 
a sonic  log  some  7 km  from  the  site;  the  log  shows  gradual  increase  in  com- 
pressional velocity  from  1.4  km/sec  at  the  surface  to  6 km/sec  at  3.6  km  depth. 

c.  Pinedale,  Wyoming;  the  site  is  located  30  km  west  of  Pinedale,  Wyoming; 
the  surface  instrumentation  consisted  of  a set  of  surface  seismographs  in  a 
vault  approximately  8 m below  the  surface  and  a microbarograph  located  15  m 
from  the  surface  vault  (Douze  and  Sherwin,  1975).  In  addition,  two  down-hole 
long-period  KS  seismograph  systems  were  installed.  The  KS  seismometer  has  not 
been  described  in  the  open  literature,  but  a brief  description  of  it  was  given 
by  Starkey  (1973).  One  was  placed  throughout  the  experiment  in  a shallow  hole 
at  a depth  of  46  m.  The  other  deep-hole  instrument  was  operated  it  different 
depths  in  an  adjacent  3000  m hole.  Because  of  an  obstruction  in  the  casing  at 
1376  m,  operation  below  this  depth  was  not  possible.  The  hole  was  filled  with 
water  below  a depth  of  190  m;  as  will  be  discussed  later,  the  presence  of 
liquid  in  the  hole  affected  the  performance  of  the  seismometer.  The  data  were 
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recorded  on  FM  analog  tape  at  sufficiently  high  levels  so  that  tape  noise  was 
not  a problem  in  the  period  range  of  20  to  100  seconds.  The  velocity  structure 
of  the  site  is  very  simple;  the  deep  hole  was  drilled  in  granite  from  a few 
feet  below  the  surface  to  the  total  depth  of  3050  m.  A sonic  and  a density 
log  were  run  and  an  almost  constant  velocity  of  5.5  km/sec  and  density  of 
2.6  g/cm^  were  measured. 


2.  THEORY 


Studies  by  Sorrells  et  al.  (1971)  and  Savino  et  al.  (1972)  have 
demonstated  that  r significant  portion  of  the  seismic  noise  with  frequencies 
lower  than  about  0.05  Hz  is  related  to  atmospheric  pressure  variations.  For 
the  purposes  of  discussion,  it  is  convenient  to  separate  the  atmospheric 
pressure  field  into  wind-generated  and  wind-free  components.  Savino  et  al . 
(1972)  were  concerned  with  the  earth  motion  in  response  to  wind-free  pressure 
changes.  In  contrast,  our  work  (Sorrells,  1971;  Sorrells  et  al.,  1971;  and 
Sorrells  and  Goforth,  1973)  has  been  concerned  primarily  with  wind-related 
earth  noise. 

In  order  to  account  for  the  observed  earth  motion,  it  is  necessary  to  make 
certain  assumptions  about  the  statistical  properties  of  the  pressure  field 
and  the  distribution  of  elastic  constants  within  the  earth.  We  originally 
assumed  that  the  wind-generated  pressure  field  could  be  approximated  by  a 
wave  which  was  convected  at  a speed  equivalent  to  the  mean  wind  speed.  This 
assumption  was  based  upon  Taylor's  Hypothesis  (Lumley  and  Panofsky,  1964) 
which  states  in  effect  that,  if  the  turbulent  velocity  fluctuations  are  small 
compared  to  the  main  stream  velocity,  the  time  variations  in  the  pressure,  as 
observed  at  a fixed  point  in  the  flow,  would  be  approximately  the  same  as  that 
due  to  the  convection  of  an  unchanging  spatial  pattern  past  the  point  with  the 
mean  flow  velocity.  In  addition,  we  initially  assumed  that  the  earth  could  be 
approximated  by  a perfectly  elastic,  homogeneous  half-space.  Later  experi- 
mental work  (Sorrells  et  al.,  1971)  demonstrated  the  need  for  a more  precise 
description  of  the  shallow  distribution  of  elastic  constants;  therefore,  a 
model  was  adopted  which  consisted  of  isotropic,  elastic  layers  separated  by 
parallel  boundaries.  A somewhat  surprising  outcome  of  that  particular  study 
was  that  once  the  layered  earth  model  was  adopted,  relatively  good  agreement 
was  obtained  between  theory  and  experiment,  despite  the  fact  that  a plane  wave 
is  a gross  oversimplification  of  the  structure  of  the  wind-generated  atmos- 
pheric pressure  field.  Sorrells  aid  Goforth  (1973)  expanded  the  theory  to 
cover  the  case  where  the  pressure  field  is  a random  process  which  is  stationary 
in  time  and  homogeneous  in  space.  In  particular,  they  derived  expressions  for 
the  transfer  function  and  the  coherence  which  relate  earth  motion  to  pressure 
changes  observed  at  a point  on  the  earth's  surface.  They  represented  the 
spatial  organization  of  the  wind-generated  pressure  field  as  scattered  waves 
traveling  at  speeds  varying  from  c = u>/(£+p)  to  c = w/(^-p)  and  arriving  from 
directions  which  vary  from  -Y  to  f,  where  t,  is  the  midpoint  and  (£±p)  the 
extrema  of  the  wave-number  range  of  interest.  For  a distributed  source  of 


this  type,  and  for  pressure  convection  velocities  small  compared  to  seismic 
velocities,  the  resultant  vertical  displacement  and  apparent  horizontal  dis- 
placement are  given  by 

V S+P 


v ^ = mr  f f Gv(k)  k dk  de 


H(6,ai)  = 


-f  S-P 


Y S+P 

kf  f • 


cos  (0-B)  Gr(k) 


k dk  de  (2) 


5-P 


where 


G (k)  is  the  vertical  component  of  the  earth's  response  to  a point 
v pressure  load 

G (k)  is  the  radial  component  (including  tilt)  of  the  earth's  response 
r to  a point  pressure  load 

k is  the  wave  number  which  is  determined  by  specifying  the 

frequency,  w,  and  the  convection  velocity,  c,  of  the  pressure 


B is  the  direction  in  which  the  tilts  are  observed 
g is  the  acceleration  due  to  gravity 


Sorrells  and  Goforth  (1973)  found  that  the  above  transfer  functions,  which 
are  based  on  an  analytical  model  possessing  all  of  the  observable  properties 
of  the  wind-generated  pressure  field,  are  experimentally  indistinguishable 
from  those  associated  with  a plane  pressure  wave.  The  analagous,  but  simpler, 
transfer  functions  for  a plane  wave  are 


V M - 

H (6,  uj)  = cos  (n-6) 
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where 


(i) 


is  the  vertical  component  of  the  earth's  response  to  a plane 
pressure  wave  of  frequency  u and  velocity  c 


0is  the  radial  component  of  the  earth's  response  to  a plane 
p - * ’ ‘ 


pressure  wave  of  frequency  id  and  velocity  c 
“ n direction  of  propagation  of  the  plane  wave 
g direction  in  which  the  tilt  is  observed 


The  plane  wave  formulations  were  used  to  compute  the  theoretical  spectra 
presented  in  this  report. 


3.  ANALYSIS  TECHNIQUES 


In  this  section  we  briefly  discuss  the  data  processing  methods  used  in  this 
study:  none  of  them  are  new  and  the  only  objective  is  to  show  how  they  were 
employed  for  the  particular  problems  associated  with  eliminating  pressure- 
generated  noise.  Appropriate  references  are  cited  for  those  interested  m 
the  details  of  any  of  the  methods. 


The  basic  inputs  required  are  auto-  and  cross-correlation  of  the  noise  recorued 
by  the  seismograph  and  microbarographs:  although  these  may  be  computed  directly, 
the  use  of  the  Fast  Fourier  Transform  (Cooley  and  Tukey,  1965)  provides  a more 
effective  approach.  The  digital  time  series  t (k)  is  divided  into  M segments 
of  K data  points;  each  segment  is  transformed  into  the  frequency  domain 


Ej  -2-nifk  \ 

yk)  exp  ( — — ) 

k=0 


(5) 


The  spectral  estimates  are  then  obtained  by  averaging  over  a number  of 
segments  (M) 

M 

PCf)  = I Y S .(f)  . S*  (f) 


(6) 


j = l 
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The  P„„  are  estimates  of  the  spectra  of  the  multiple  time  series  and  Pmm  are 
the  cross  spectra.  The  asterisk  indicates  complex  conjugate.  This  approach 
to  the  computation  of  power  spectra  was  discussed  in  detail  by  Welch  (196"). 

At  this  point  it  is  possible  (Burg,  1964)  to  compute  optimum  filters  directly 
in  the  frequency  domain.  For  reasons  that  will  be  apparent  later,  we  prefer 
to  use  the  inverse  Fourier  Transform  in  the  results  of  equation  (6)  and  thus 
obtain  the  auto-  and  cross-correlations  of  the  time  series. 


The  optirii'im  filters  used  have  been  described  (see  Robinson,  1967);  we  will 
only  briefly  summarize  their  applicability  for  the  purpose  of  this  study. 
Suppose  we  have  N input  microbarographs  y^  (t),  i = 1,  N and  we  wish  to 
predict  the  amount  of  energy  on  an  output  channel,  x(t),  the  earth  noise  on 
the  seismograph  that  is  caused  by  the  pressure  variations  recorded  by  the 
microbarographs.  We  desire  a set  of  N filters,  f(t)  of  length  T,  that  can 
be  applied  to  the  N microbarographs  to  obtain  an  output  z(t)  that  will  closely 
approximate  the  desired  output  channel  x(t). 


N L 

z(tj  = £ £ fA(s)  • yi(t-s) 
i=l  s=0 


s = 0,  L 
t = 0,  T + L + 1 


(7) 


We  seek  the  filters  by  minimizing  the  error  function  defined  as 


E 


T+L+l 


£ 


_ 2 

N L 

£ £ fj(s)  yi(t-s)  - x(t) 
i=l  s=0 


(8) 


The  error  E is  minimized  when  the  partial  derivatives  with  respect  to  the 
filter  coefficients  are  put  equal  to  zero 


BE/8fk  (s)  = 0 for  k = 1 , N;  s = 0,1. 

The  derivation  is  straightforward  (Robinson,  1967)  and  results  in: 

N L 

^ 2 Rki  (u"s)  fi(s)  = ck(u)  O) 

i=l  s=0 

where  R^  is  the  multichannel  correlation  matrix  of  the  inputs,  the  microbaro- 
graph channels  we  are  using,  and  the  Cj<  are  the  cross-correlations  between  the 
desired  output,  the  seismograph,  and  the  microbarographs. 

Examination  of  the  equation  (9)  shows  that  it  consists  of  a large  number  of 
simultaneous  equations.  Fortunately,  there  is  available  an  effective  recursive 
method  of  solving  for  the  desired  filters;  e.g.,  we  start  by  computing  N 1- 
point  filters  and  from  them  compute  a 2-point  filter,  and  so  forth.  The 
recursive  approach,  first  described  by  Levinson  (1949),  is  made  possible  by 
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the  fact  that  the  autocorrelation  matrix  is  a Toeplitz  matrix.  The  advantage 
of  the  approach  is  that  at  each  reiteration  the  remaining  error  is  computed 
and  the  effectiveness  of  that  length  of  filter  is  known.  Uusally  the  error 
will  decrease  rapidly  and  then  level  off,  thus  one  finds  the  least  number  of 
filter  points  required.  When  processing  large  amounts  of  data,  the  savings 
in  processing  can  become  appreciable. 

When  processing  data  off-line,  one  is  not  confined  to  realizable  filters 
(i.e.,  using  only  past  values)  but  can  also  use  future  data  points.  In  some 
cases  the  capability  of  the  optimum  filters  is  greatly  improved  by  using  non- 
realizable  filters.  However,  because  of  the  phase  delays  involved  between 
the  seismograph  ard  microbarograph  traces,  no  improvement  is  obtained  by  this 
procedure  for  the  case  under  study  here. 


By  means  of  coherence  measurements  it  is  possible  to  obtain  an  estimate  of  the 
fraction  of  the  power  that  can  be  predicted.  The  ordinary  coherence  function 
(or  coherence  scared)  is  obtained  from  the  auto-  and  cross-power  spectra 
in  the  usual  fashion 


I ^mn  (O  1 2 

^mm(f)  ^nn^O 


(10) 


It  is  a measure,  frequency  by  frequency,  of  how  much  of  the  power  on  a 
seismograph  trace  can  be  predicted  from  the  microbarograph  input. 


In  the  same  fashion,  the  multiple  coherence  (Bcndat  and  Picrsol,  1966)  can  be 
used  for  the  multichannel  case;  it  is  defined  from  the  spectral  matrix  <t>ff) 
of  all  the  traces 


J *jj(f)  *jj(f) 

4>jj(f)  = power  spectrum  of  the  jth  trace,  the  seismograph 

(f)  = diagonal  element  of  the  inverse  of  the  spectra)  matrix. 


The  multiple  coherence  is  a measure  of  the  linear  relationship  between  an 
output  in  this  case  seismograph,  and  a number  of  inputs,  the  microbarographs. 
As  such,  it  is  a measure  of  how  well  a multichannel  filter  would  be  able  to 
predict  the  noise  recorded  by  the  seismograph.  As  it  is  more  economical  to 
compute,  we  have  used  it  extensively  instead  of  actually  computing  the  multi- 
channel filter  and  applying  it  to  the  data.  One  problem  with  the  multiple 
coherence  is  the  rapid  increase  of  the  bias  in  the  estimate  with  increasing 
number  of  inputs  for  low  coherences  (White,  1973). 
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Extensive  use  was  made  of  frequency -wave  number  spectra  in  order  to  study  the 

structure  of  atmospheric  pressure  variations.  Our  method  for  making  the!® 

estimates  is  an  adaptation  of  an  efficient  computational  algorithm  described 

by  Smart  and  Flinn  (1971).  The  basic  change  was  to  extend  the  block  averagi  g 

tprhninne  (Welch  1967)  to  the  estimation  of  frequency  wave -number  spec  ra. 
technique  (Welch , iw/j  to  uie  Vnttr.ilkr  transform  of  the  £th  time 


technique  iweicn,  i»o/j  tu  me  --  --  - _ , - 

Briefly,  suppose  that  F-j  8,(01)  is  the  finite  Fourier  transform  of  the  l 

segment  recorded  by  the  jth  detector.  It  may  be  written  as 


Fjj^w)  = Ajt(u)  { exp l - ioc j p (u>)  ] ) 


(12) 


The  power  spectral  and  cross-power  spectral  density  estimates  are  then  given 
by 

V“>  = t E ,1') 

t-l 


where  * denotes  complex  conjugation.  Equation  (13)  yields  power  spectral 
density  estimates  for  m = j and  cross-power  spectral  density  estimates  for 

m t j . 

The  frequency  wave-number  spectral  estimate  is  defined  as 


V V 

P(to,k)  = E E Sjm(w)  explik  • (rj  - rm)  1 

j = l m=l 


(14) 


• . . ,1:  ot nnre  to  the  i th  detector  with  reference  to  an 

SHraiy'orlgln  ^e„  Jhrough  the  use' of  equations  (13)  and  (.3,  (Sorrells 
and  Douze,  1974) 

' 2 


P (ui , k ) = -r 

L £.  = 1 


N 


52  Ajt(u)  exp(-iajfc(u)  ♦ ik  * rj 

j = l 


(15) 


Basically  (15)  indicates  that  the  f-k  estimates  for  each  time  segment  can  be 
added . The  Smart  and  Flinn  algorithm  is  used  to  obtain  these  estimates  and 
the  Fisher  statistic  is  computed  (Smart  and  Flinn,  1971 J in  order  to  as.-  g 
statistical  significance  to  the  results. 


\ 
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4.  SOURCES  OF  SEISMIC  NOISE  IN  THE  20-100  SECOND  PERIOD  RANGE 


The  noise  recorded  by  a long-period  seismograph  can  be  broadly  classified  as 
either  earth  noise  or  "system"  noise.  Earth  noise  is  the  result  of  ambient 
vibrations  of  the  ground  regardless  of  whether  the  cause  is  elastic  waves  or 
quasi-static  deformations  caused  by  pressure  variations  in  the  atmosphere. 
"System"  noise  is  used  in  this  report  both  for  the  noise  generated  in  the 
seismograph  system  and  the  noise  within  the  immediate  environment,  principally 
the  vault.  The  noise  from  the  immediate  environment  is  difficult  to  separate 
from  other  sources:  effects  such  as  buoyancy  (Crary  and  Ewing,  1952),  convec- 

tion currents  in  casing  or  vault  or  bending  of  base  plates  (Holcomb,  1975), 
can  all  contribute  to  the  environmental  noise  and  are  especially  severe  on 
surface  horizontal  seismographs. 

The  system  noise  caused  by  the  seismograph  system  is  better  known.  Magnifi- 
cation at  which  long-period  seismographs  are  presently  operated  is  such 
that  the  noise  from  thermal  agitation  approaches  the  earth  noise.  A number 
of  authors  have  investigated  theoretical  and  experimental  noise  levels. 

Almost  all  experimental  data  are  for  standard  mass-spring  systems;  the  system 
noise  of  a force-balance  seismograph  (KD  seismometer)  is  more  difficult  to 
determine  because  the  usual  method  of  replacing  the  mass  with  an  equivalent 
resistor  is  not  applicable. 

Theoretical  system  noise  levels  have  been  recently  discussed  in  detail  by 
Fix  (1972)  and  Melton  (1975).  Experimental  results  have  been  shown  by  Fix 
(1972),  and  Savino  (1971)  for  conventional  seismographs  and  by  Douce  and  Starkey 
(1973)  for  the  force-balance  system  used  at  Pinedale.  In  all  cases  the  system 
r.oise  is  sufficiently  low  not  to  seriously  affect  the  experimental  results; 
however,  it  represents  an  appreciable  percentage  of  the  total  noise  recorded. 
Furthermore,  it  should  be  noted  that  system  noise  levels  are  almost  always 
measured  after  careful  maintenance  of  the  system,  and  that  the  system  ncise  may 
be  appreciably  higher  under  more  normal  operating  conditions. 

The  background  noise  at  periods  greater  than  20  seconds  is  quite  stable  as 
compared  to  the  microseismic  peaks  which  vary  up  to  two  orders  of  magnitude. 
Spectral  estimates  of  the  noise  recorded  in  the  20-100  second  period  range 
by  the  vertical  seismographs  are  shown  in  figure  1.  The  data  from  which  these 
estimates  were  made  were  taken  during  calm  intervals  in  the  month  of  January 
1974,  at  McKinney,  Texas.  Observe  that  the  spread  of  the  estimates  is  about 
3 dB.  The  anticipated  spread  of  estimates  at  the  90  percent  confidence  level 
for  a time  stationary  process  is  2.8  dB.  Thus,  these  results  imply  that  in 
the  absence  of  wind-generated  earth  motion,  the  vertical  seismic  noise  may  be 
considered  to  be  time  stationary  for  intervals  of  time  at  least  on  the  order 
of  a month.  These  results  are  consistent  with  observations  reported  by 
Savino  et  al.  (1972)  regarding  the  spectra  of  seismic  noise  observed  in  the 
Ogdensburg  mine  at  a depth  of  542  meters.  Agreement  in  this  area  is  to  be 
expected  since  the  spectra  observed  at  depths  in  excess  of  several  hundred 
meters  should  be  free  of  wind  related  earth  noise  most  of  the  time.  There  is 
also  some  evidence  which  suggests  that  the  calm  interval  vertical  seismic  noise 
spectrum  may  vary  slightly  from  winter  to  summer  at  McKinney,  Texas,  as  seen 
when  the  mean  of  the  January  1974  estimates  is  compared  to  the  estimates  of 
September  1973. 
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ESTIMATED  POWER  SPECTRAL  DENSITY  (mM2/Hz) 


PERIOD  (seconds) 

Figure  1.  A comparison  of  noise  power  spectral  densities.  Data  were 

recorded  by  a vertical  seismograph  during  intervals  when  the 
mean  wind  speed  was  less  than  1 m/sec.  Spectral  estimates 
have  been  corrected  for  system  response,  McKinney,  Texas. 


G 7569 
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The  September  1973  estimate  is  about  3 dB  lower  than  the  January  mean.  This 
result  suggests  that  the  10  dB  seasonal  variation  in  spectral  amplitudes 
observed  by  Savino  et  al.  (1972)  is  strongly  suppressed  at  the  McKinney  site. 

On  the  other  hand,  the  results  at  Pinedale,  Wyoming  varied  approximately  10  dB 
over  any  month  during  the  course  of  the  experiment,  June  through  December  .974. 

Considering  the  numerous  noise  sources,  i.e.,  system  and  environmental,  it  is 
often  difficult  to  determine  whar  is  truly  "earth  noise."  If  we  define  earth 
noise  as  that  portion  of  the  noise  that  is  coherent  over  distances  of  the 
dimension  of  the  vault,  an  estimate  of  earth  noise  can  be  obtained  by  operating 
tvo  seismographs  in  close  proximity.  The  ratio  of  earth  noise  to  other  noise 
sources  can  be  calculated  from  the  power  spectral  density  estimate  and  the 
coherence  estimates  between  the  outputs  of  the  two  seismographs.  Under  the 
appropriate  experimental  conditions  the  coherent  power  can  be  attributed  to 
earth  noise  and  the  incoherent  power  to  system  noise.  Let  and  p7  denote 
the  ESR  (earth/system  noise  ratio)  for  each  of  the  two  seismographs  and  let 
Y^2  represent  the  square  of  the  coherence  between  their  outputs.  Then  it  can 
be  shown  that , 


y<=  = 

12  1 ♦ Pi 

(e.g.  Foster  and  Guinzy,  1967) 


p2 

1 ♦ Pj 


(16) 


The  ratio  of  power  spectra  R2  can  also  be  written  in  terms  of  p1  and  p 
being  given  by,  12  Z 


R2 

12 


1 ♦ Pi 


1 ♦ p- 


Equations  (16)  and  (17)  may  be  solved  for 


O 

' 1 


and  p^  yielding 


(17) 


Y12 


P,  = 


1 U12  _ Y12 


P2  = 


'12 


12 


12 


(18) 


For  the  surface  vertical  seismographs  at  McKinney  it  has  been  found  that 
Pi  « P2*  Typical  values  estimated  from  calm  interval  samples  are  shown  in 
figure  2.  For  reference,  we  have  also  included  values  of  the  ratio  estimated 
for  the  vertical  seismograph  in  the  mine  at  Grand  Saline,  Texas,  (Sorrells 
et  al.,  1971'.  It  will  be  observed  that  at  periods  greater  than  about 
FO  seconds  tne  ratios  at  McKinney  are  lower  than  those  at  Grand  Saline.  The 
cause  of  this  discrepancy  is  not  clearly  understood  at  the  present  time.  It 
appears  to  be  related  to  a difference  in  the  amplifiers  used  at  the  two  loca- 
tions. A phototube  amplifier  with  a 30-second  galvanometer  was  used  at  Grand 
Saline.  A solid-state  amplifier,  instead  of  a phototube  amplifier,  is  being 
used  at  McKinney.  The  lower  ESR  values  found  during  our  current  program  do 
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not  seriously  influence  the  outcome  of  our  investigation.  They  do,  however, 
limit  the  resolution  with  which  pressure  related  earth  motion  can  be  separated 
from  the  noise  generated  by  other  sources,  particularly  at  periods  greater  than 
60  seconds. 

During  windless  periods  at  McKinney,  fr  quency-wave  number  spectra  often 
show  significant  peaks  (from  the  Fishe.-  statistic)  at  high  velocities.  Because 
the  microbarograpn  array  used  is  small,  it  is  difficult  to  determine  whether 
the  peaks  are  associated  only  with  infrasonic  waves  traveling  at  330  m/sec. 

An  example  of  significant  energy  traveling  at  infrasonic  velo:ities  is  shown 
in  figure  3.  However,  a significant  number  of  f-k  spectra  show  even  higher 
velocities  such  as  shown  in  figure  4.  The  origin  of  these  high-velocity  pres- 
sure variations  is  not  understood  but  they  occur  at  levels  sufficiently 
high  vo  be  recorded  by  the  seismographs.  The  RMS  amplitude  range  of  these 
pressure  changes  is  generally  0.1  to  5 pbars.  At  soft  rock  sites  in  the  period 
range  pressure  variations  of  this  order  can  be  caused  by  winds  with  mean 
velocities  as  low  as  0.5  m/sec. 

In  the  f-k  analysis  the  presence  of  wind  noise  rapidly  decorrelating  with 
distance  has  the  same  effect  as  adding  white  noise  to  the  spectrum  and  the 
coherence  between  a single  microbarograph  and  a seismograph  will  be  significantly 
decreased.  The  reason  for  this  is  that  the  square  of  the  ordinary  coherence  is 
a measure  of  the  seismic  noise  power  that  is  linearly  predictable  from  a single 
microbarograph.  However,  as  a general  rule,  not  all  of  the  pressure  related 
earth  noise  recorded  by  a seismograph  is  predictable  from  the  output  of  a co- 
located microbarograph.  The  sole  exception  occurs  when  variations  in  the 
atmosphereic  pressure  are  the  result  of  a plane  wave  (Sorrells  and  Goforth, 

1973).  As  the  structure  of  the  pressure  variations  depart  from  this  simple 
form  the  percentage  of  predictable  earth  noise  will  decline  ever  though  the 
total  pressure  contribution  may  remain  constant.  For  this  reason,  it  is 
sometimes  necessary  to  resort  to  multiple  coherence  estimates  to  detect  the 
presence  of  a pressure  related  component.  In  our  case,  the  square  of  the 
multiple  coherence  is  a measure  of  the  percentage  of  the  seismic  noise  power 
that  is  predictable  from  the  outputs  of  an  array  of  microbarographs.  In  the 
case  of  poorly  organized  pressure  fields,  the  percentage  of  power  predictable 
from  a microbarograph  array  should  be  a better  approximation  to  the  total 
pressure  related  seismic  noise  power  than  that  obtained  from  the  ordinary 
coherence.  However,  in  the  case  of  well  organized  fields,  ordinary  and  multi- 
ple coherence  estimates  should  coincide  ar.d  should  be  equal  to  the  percentage 
of  seismic  noise  power  that  is  caused  by  atmospheric  pressure  variations. 

One  of  the  more  interesting  results  of  our  investigations  to  date  is  that 
close  agreement  between  ordinary  and  multiple  coherence  estimates  is  not  un- 
common. Examples  calculated  from  data  recorded  during  a calm  interval  on 
29  January  1974,  are  shown  in  figure  5.  Frequency  wave-number  estimates  of 
variations  in  the  atmospheric  pressure  field  for  the  same  interval  yielded 
somewhat  ambiguous  results  because  of  the  relatively  small  aperture  of  the 
microb  rograph  array.  There  is,  however,  some  suggestion  that  infrasonic 
waves  ace  the  principal  source  of  atmospheric  pressure  variations  during  this 
interval.  (See  figure  3.)  Suppose  we  assume  that  variations  in  the  atmos- 
pheric pressure  field  during  this  interval  indeed  are  the  result  of  scattered 
infrasonic  waves  whose  speed  range  is  bracketed  by  the  values  and  cu.  If 
the  field  is  stationary  in  time  and  homogeneous  in  a plane  parallel  to  the 
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Distribution  of  pressure  on  the  wave  number  plane  at 
period  of  26.95  seconds.  Results  indicate  a source 
infrasonic  atmospheric  pressure  variations  from  the 
northwest,  McKinney,  Texas. 
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earth's  surface  and  the  power  spectra  of  the  waves  are  independent  of  their 
speeds,  then  for  the  case  of  a vertical  seismograph  and  a microbarograph 
located  at  the  surface  of  a homogeneous  and  isotropic  half  space,  the  square 
of  the  ordinary  coherence  is  given  by 


(“) 


$(oj)  (19) 

(Sorrells  and  Goforth, 

1973,  eq.  61) 


where  $ (co)  is  the  percentage  of  seismic  noise  power  caused  by  atmospheric 
pressure  variations.  The  term  ii  brackets  is  the  ratio  of  predictable  pressure 
related  noise  power  to  total  pressure  related  noise  power. 


This  ratio  decreases  monotonically  as  the  cu  increases  with  respect  to  cz  (i.e., 
as  the  speed  range  of  the  scattered  waves  expands) . Now  one  would  expect  the 
speed  range  of  scattered  infrasonic  waves  to  be  no  greater  than  about  300-600 
meters/second  / cu  < 2\. 

) 

For  a speed  range  of  this  order  the  ratio  of  predictable  to  total  power  is 
greater  than  0.9.  Thus,  from  equation  (19)  the  square  of  the  ordinary  co- 
herence will  underestimate  the  percentage  of  pressure  related  noise  power  by 
10  percent  or  less.  Since  the  square  of  the  multiple  coherence  is  always 
greater  than  or  equal  to  the  square  of  the  ordinary  coherence  but  less  than 
or  equal  to  the  percentage  of  pressure  related  noise  power  it  follows  that  in 
the  case  of  scattered  infrasonic  waves  the  ordinary  and  multiple  coherences 
should  be  approximately  the  same.  It  is  therefore  our  belief  that  the  results 
in  figure  5 together  with  data  from  the  frequency-wave  number  estimates  indi- 
cate that  infrasonic  waves  are  the  principal  source  of  the  atmospheric  pressure 
variations  and  pressure  related  earth  noise  observed  during  this  interval. 


It  should  be  noted  that  relatively  very  minor  amounts  of  wind- generated  noise 
are  sufficient  to  obscure  the  infrasound-associated  relationship  between  the 
seismographs  and  microbarographs.  But  in  the  seismic  background  spectra  the 
percentage  of  pressure-related  noise  caused  by  wind  is  less  than  10  percent 
of  the  total  of  such  noise,  except  when  the  power  of  the  wind  pressures  exceeds 
that  of  the  infrasonic  component  by  fully  a factor  of  20  or  more  in  the 
microbarograph  spectra.  This  illustrates  the  point  that  though  the  contribution 
of  infrasonic  waves  to  the  pressure  spectra  may  be  relatively  small,  it  may  still 
account  for  virtually  all  of  the  pressure  related  earth  noise  found  in  the 
vertical  seismic  noise  spectra. 


When  examining  the  numerous  spectra  of  the  noise  obtained  from  the  seismometers 
at  Tinedale,  where  the  seismographs  are  unaffected  by  the  wind,  it  is  clear 
that  at  times  there  is  considerable  structure  in  the  noise  in  the 
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notch  observed  at  20-40  seconds  of  period.  Small  but  statistically  significant 
peaks  are  present;  the  reasons  for  believing  they  consist  of  seismic  energy  are: 

a.  They  appear  on  both  horizontals  and  verticals  at  about  the  same 
power  in  the  same  fashion  as  the  microseismic  peaks  at  18-  and  6-second 
periods.  If  they  were  associated  with  environmental  noise  one  would  expect  them 
to  be  considerably  larger  on  the  horizontals,  which,  because  of  their  tilt 
sensitivity,  almost  invariably  are  more  susceptible  to  environmental  noise. 

b.  Associated  with  the  peak  in  the  spectra  there  is  usually  an  increase 
in  coherence  between  instruments  of  ' n same  orientation  operating  at  different 
depths . 

c.  They  do  not  appear  to  be  related  to  the  atmosphere  as  there  is  never 
a corresponding  peak  in  the  power  spectra  of  the  microbarograph  or  any  increase 
in  the  coherence  between  any  seismometer  and  the  microbarograph.  Furthermore, 
when  peaks  in  quiet  day  spectra  were  found  at  McKinney,  the  frequency-wave 
number  spectra  did  not  show  any  long-wave  length  atmospheric  phenomena  at  the 
corresponding  periods. 

It  will  require  an  array  of  high-quality,  burned,  long-period  seismographs 
such  as  is  presently  being  constructed  in  Iran  to  verify  the  conclusion  that 
there  is  appreciable  seismic  energy  at  periods  of  20  to  60  seconds. 


5.  PREDICTION 


The  experiments  in  predicting  long-period  noise  were  carried  out  at  McKinney, 
Texas,  because  the  required  microbarographs  were  available  at  this  site.  Un- 
fortunately, even  during  windless  intervals,  only  the  vertical  long-period 
seismographs  could  be  operated  at  magnifications  comparable  to  those  attainable 
at  some  depth  below  the  surface.  The  horizontal  seismographs  are  noisy  even 
during  windless  periods  and  the  prediction  filtering  techniques  for  the  hori- 
zontals were  far  less  effective  than  for  the  vertical  seismographs. 


5.1  WINDY  PERIODS:  VERTICAL  SEISMOGRAPHS 

Amplitudes  of  wind- induced  pressure  variations  can  easily  be  1 to  2 orders  of 
magnitude  greater  than  long  wavelength  infrasonic  or  other  atmospheric  phenomena 
detected  during  windless  peiiods.  Therefore,  during  intervals  of  moderately 
high  wind  speeds  (>  4 m/sec),  the  pressure-related  noise  will  be  dominated  by 
quasi-static  deformations  in  response  to  the  wind-generated  pressure  variations. 
Because  of  the  short  wavelengths  associated  with  wind-induced  turbulence, 
correlation  is  found  only  between  the  outputs  of  the  seismographs  and  -he  co- 
located microbarograph.  Tt  is  therefore  unnecessary  to  include  the  outputs  of 
the  other  microbarographs  in  the  processing;  the  experimental  problem  thus 
simplifies  to  the  estimation  of  single-channel  optimum  filters. 
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In  the  following,  note  that  none  of  the  spectra  of  this  section  are  corrected 
for  system  response. 

Figure  6 shows  the  results  of  one  of  the  experiments;  shown  are  the  spectra 
of  the  vertical  earth  noise  as  recorded  by  the  seismograph  and  the  noise  after 
the  predictable  portion  of  the  noise  had  been  subtracted.  The  wind  was  blowing 
approximately  8 m/sec  during  this  time.  The  required  auto-  and  cross- 
correlation were  computed  from  a 2-hour  sample  and  the  calculated  filter  applied 
to  the  same  time  interval.  The  filter  used  was  58  points  (232  seconds)  long 
and  had  no  delay  built  in;  several  filters  with  different  delays  were  tried 
but  no  improvement  could  be  obtained  by  using  a non-realizable  filter.  Exami- 
nation of  figure  6 shows  clearly  that  a large  percentage  of  the  earth  noise 
is  caused  by  wind-induced  pressure  variations  because  this  noise  was  effective- 
ly predicted  by  an  optimum  filter  based  on  the  microbarograph  output.  In  order 
to  show  the  effectiveness  of  the  process,  the  earth  noise  on  a windless  day 
(about  a week  before)  is  also  included.  It  has  been  found  that  in  the  absence 
of  wind  related  earth  noise,  the  spectrum  of  the  vertical  component  of  the 
noise  remains  approximately  stationary  (within  6 dB)  foi  intervals  at  least  on 
the  order  of  a month.  As  expected,  at  periods  of  20  seconds  and  less  the 
filter  is  ineffective  because  the  microseisms  predominate  in  the  earth  noise. 
Examination  of  tne  mean  square  error  shows  that  55  percent  of  the  total  power 
was  predictable  with  this  filter.  The  total  power  includes  the  microseisms; 
in  the  band  from  20-  to  100-second  periods,  the  percentage  of  predicted  energy 
is  considerably  greater. 

Figure  7 shows  a plot  of  average  wind  velocity  versus  the  percentage  of  total 
energy  that  can  be  predicted  with  a single-channel  filter.  The  filters  were 
all  computed  from  1-hour  sample  times  when  the  wind  was  fairly  constant.  While 
there  is  considerable  scatter  in  the  data  points  there  is,  as  expected,  a 
distinct  increase  of  prediction  capability  with  increasing  level  of  wind  and 
therefore,  atmospherically- generated  earth  noise.  The  scatter  of  the  points 
is  at  least  partly  caused  by  variations  of  the  power  level  of  the  microseisms. 
Despite  the  variability  shown  in  figure  7 the  stability  of  the  computed  filters 
is  good,  as  will  be  discussed  later  in  this  section. 

The  single  channel  optimum  filters  discussed  in  the  paragraphs  above  were 
estimated  from  data  recorded  during  a finite  interval  of  time.  They  were 
then  used  with  considerable  success  to  reduce  seismic  noise  levels  over  the 
same  interval.  However,  in  practice,  it  is  certainly  desirable  and  some- 
times necessary  to  apply  the  filters  computed  for  one  interval  to  the  data 
recorded  during  other  intervals  as  well.  Information  regarding  the  temporal 
stability  of  the  filter  estimates  is,  therefore,  of  considerable  practical  im- 
portance. In  order  to  obtain  such  information  we  computed  filters  from 
seventy-six  1-hour  samples  spanning  a 2-week  interval.  Typical  estimates  of 
the  modulus  and  phase  of  the  filter  transfer  functions  are  shown  in  figures 
8 and  9.  These  results  were  obtained  from  data  recorded  during  oix  of  the 
intervals  when  the  mean  wind  velocity  was  6-7  m/sec  from  the  south.  The 
dashed  curves  bracket  the  90  percent  confidence  interval  for  sample  78-0500. 
These  results  imply  that  the  observed  scatter  can  be  attributed  to  normal 
statistical  variability.  It  is  important  to  note  that  the  data  from  which 
these  filters  were  estimated  were  taken  only  during  intervals  when  the  mean 
wind  speed  and  direction  were  approximately  the  same.  Theoretical  considera- 
tions (Sorrells,  1971),  lead  us  to  expect  systematic  changes  in  the  filters 
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POWER  SPECTRUM  (m/j 2/Hz)  RELATIVE  TO  40  SECOND  PERIOD 


Figure  6.  Power  spectra  of  the  noise  recorded  by  a vertical  seismograph 

(a)  original  recording,  (b)  after  signal  channel  optimum  filter, 
(c)  quiet  day  spectrum.  Wind  velocity  8 m/sec,  McKinney,  Texas. 
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Figure  7.  Average  wind  velocity  versus  the  percentage  of  total  power 
predicted  by  a single-channel  filter,  McKinney,  Texas. 
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Moduli  of  the  transfer  functions  of  the  optimum  filters 
used  to  reduce  noise  on  the  vertical  seismograms, 
McKinney,  Texas. 
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as  the  mean  wind  speed  and  direction  changes.  For  example,  in  the  case  of  the 
homogeneous  and  isotropic  half  space,  the  modulus  of  the  vertical  filter  should 
increase  in  proportion  to  the  mean  wind  speed  while  the  phase  remains  constant. 
On  the  other  hand,  the  modulus  of  the  horizontal  filter  should  be  independent 
of  the  wind  speed  but  should  vary  as  the  cosine  of  the  angle  between  the 
sensitive  axis  of  the  seismograph  and  wind  direction.  In  aalition  the  phase 
should  remain  constant  as  long  as  the  wind  direction  remain,  within  + 90  degrees 
of  the  positive  direction  of  the  sensitive  axis.  If  it  exceeds  these  bimits, 
the  phcse  of  the  filters  should  shift  by  180  degrees. 

In  order  to  detect  the  anticipated  changes  in  the  vertical  filter,  we  have 
computed  the  average  modulus  and  phase  for  the  4-6,  6-8,  and  8-10  m/sec  wind 
speed  ranges.  These  are  shown  in  figures  10  and  11.  As  expected  the  modulus 
increases  as  the  wind  speed  increases  while  the  phase  remains  constant.  Our 
data  base  is  not  currently  adequate  to  demonstrate  the  cosine  dependence  of 
the  modulus  of  the  horizontal  filter,  as  the  wind  was  either  out  of  the  north 
or  the  south  during  the  2-week  interval  selected  for  detailed  study.  However, 
the  preliminary  results  suggest  that  the  modulus  of  the  horizontal  filter  is 
independent  of  the  wind  speed  as  implied  by  the  theory. 

On  the  basis  of  these  results,  we  may  infer  that  the  optimum  filters  will  re- 
main stable  as  long  as  the  mean  wind  speed  and  direction  remain  approximately 
stationary.  It  has  been  our  experience  that,  apart  from  the  situation  during 
frontal  passages,  these  quantities  tend  to  vary  slowly  and  smoothly  with  time. 

In  the  north  Texas  area,  for  example,  the  wind  direction  may  remain  roughly 
constant  for  periods  of  time  on  the  order  of  days,  while  the  mean  wind  speed 
undergoes  a smooth  diurnal  variation.  It  is  also  important  to  observe  that 
the  filters  for  the  vertical  seismograph  vary  in  a predictable  fashion  with 
respect  to  the  mean  wind  speed  and  direction.  These  results  suggested  that 
time  varying  filters  whose  properties  change  in  response  to  changes  in  the 
mean  wind  could  be  used  to  effectively  remove  wind-generated  noise.  However, 
a few  preliminary  tests  of  this  approach,  using  adaptive  filters  based  on 
the  Widrow  algorith  (Widrow,  1968),  were  not  successful,  primarily  because 
of  instability  in  the  process.  Further  investigations  are  required  to  show 
whether  such  an  adaptive  filtering  technique  based  directly  on  mean  wind 
speed  and  direction  will  prove  superior  to  the  one  used  above. 

To  present  examples  of  the  improvements  obtained  on  seismograms,  visually, 
which  is  the  way  most  seismic  traces  are  analyzed,  we  show  (figure  12)  a 
couple  of  representative  cases,  one  for  noise  only  on  a horizontal  and  one  for 
a vertical  with  a signal.  In  each  of  the  two  cases,  the  top  trace  is  the 
original  recording.  The  middle  trace  shows  the  output  of  a single-channel 
prediction  filter,  and  the  bottom  trace  is  the  result  of  subtracting  the  pre- 
dicted time  series  from  the  original  recording.  In  both  cases,  the  wind  was 
blowing  8-10  m/sec,  and  the  auto-  and  cross-correlations  for  the  filter  design 
were  estimated  from  1 hour  of  data.  The  noticeable  phase  shift  between  the 
original  and  predicted  traces  is  a function  of  the  particular  data.  In  other 
parts  of  the  same  records  these  slight  phase  shifts  are  different.  Figure  12b 
shows  the  results  of  applying  the  technique  to  the  surface  wave  group  associa- 
ted with  an  earthquake  with  an  m^  of  3.6  that  occurred  south  of  Panama 
approximately  31  degrees  south  of  our  site.  Notice  how  the  onset  of  the 
surface  waves  on  the  observed  vertical  seismograph  has  been  obscured  by  the 
almost  simultaneous  arrival  of  a wind-generated  noise  pulse  shown  on  the 
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Figure  10.  comparison  of  the  average  moduli  of  the  vertical  optimum 
filters  for  various  mean  wind  speeds,  McKinney,  Texas. 
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Figure  12.  Observed  and  optimum  filtered  seismograms  (a)  noise 
recorded  by  horizontal  seismograph  (LPN) , (b)  signal 
recorded  by  vertical  seismograph  (LPZ) . 
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predicted  trace.  Subtraction  of  the  predicted  trace  from  the  observed 
seismogram  works  very  well  in  cleaning  up  the  early  cycle  or  two  of  the  sur 

face  wave  group. 


5.2  CALM  PERIODS:  VERTICAL  SEISMOGRAPHS 

During  windy  periods,  single-channel  optimum  filtering  of  the  vault 
microbarogram  is  effective  in  reducing  its  noise.  During  windless  periods 
(<1.0  m/sec)  the  nature  of  the  pressure  field  is  far  more  complex.  As  dis- 
cussed in  section  4,  frequency-wave  number  analysis  of  the  microbarograph 
array  shows  the  presence  of  multiple  sources  and  a wide  range  of  velocities. 

The  array  is  not  sufficiently  large  to  clearly  differentiate  between  ti.ese 
sources  and  velocities,  for  it  was  designed  to  study  the  6- second  microbaroms. 
At  present,  it  is  not  possible  for  us  to  tfake  definite  statements  about  the 
nature  of  the  pressure  field  on  quiet  days. 

Because  computation  of  multichannel  filters  is  time  consuming,  the  multiple 
coherence  in  a given  case  was  usually  calculated  to  determine  beforehand 
if  the  optimum  filters  would  be  effective.  The  results  varied  considerably. 
Figure  13  shows  two  representative  examples.  As  mentioned  in  the  section  on 
ata  processing,  the  bias  problem  of  multiple  coherence  is  severe,  shown  on 
the  figure  is  the  expected  zero  coherence  level  computed  accordin'’  to  the 
maximum  likelihood  technique  (White,  1973).  The  samples  were  tak ..  5 days 
apart  and  only  six  micro>>arographs  were  used  because  of  instrumentation  prob- 
lems with  number  2. 

It  is  apparent  from  these  results  (figure  13)  that  there  would  be  a large 
difference  in  effectiveness  of  the  multichannel  filters  for  these  two 
intervals.  For  1 day  the  multiple  coherence  is  essentially  negligible  except 
for  one  peak  at  39  seconds,  and  multichannel  filters  would  not  predict  the 
seismograph  noise.  For  the  other  day  the  multiple  coherence  is  significantly 
above  the  expected  zero  level  for  the  periods  between  40  and  128  seconds;  thus 
the  filters  would  effectively  reduce  the  noise  level.  This  second  curve 
closely  resembles  the  multiple  coherences  of  the  data  used  to  compute  the 
optimum  filters  described  below. 

Figure  14  shows  an  example  of  optimum  filtering  during  a quiet  day  when  the 
coherence  measurements  indicated  clearly  that  multichannel  filtering  would  be 
effective,  while  the  coherence  between  the  seismograph  and  the  vault  micro- 
barograph showed  that  a single-channel  filter  reduced  the  power  level  only 
marginally.  The  results  are  quite  clear;  however,  the  improvement  is  not 
as  great  on  this  quiet  day  as  it  is  for  intervals  with  wind-generated  noise. 

The  sources  in  the  atmosphere  that  were  contributing  to  the  pressure- induced 
noise  on  the  seismograph  at  the  time  these  filters  were  computed  were  proba- 
bly low  level  infrasonic  waves.  The  frequency-wave  number  spectrum  shows 
that  the  energy  between  50-60  second  periods  is  at  least  partially  the  result 
of  acoustic  waves  from  the  northwest.  However,  the  large  difference  in  pre- 
dictability from  single-  to  four-channel  filters  clearly  indicates  that  the 
acoustic  waves  are  not  the  only  source  of  pressure  variations. 
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jnples  of  multiple  coherences  obtained  on  windless 
Vertical  seismograph  at  McKinney,  Texas. 
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Power  spectra  of  the  noise  on  a vertical  seismograph  during 
a windless  day.  (a)  original  recording,  (b)  single-channel 
filter,  (c)  4-channel  optimum  filters,  McKinney,  Texas. 
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Out  of  the  total  power  recorded  by  the  seismograph,  the  single-channel  filter 
could  only  predict  13  percent;  the  four-channel  set  of  optimum  filters  pre- 
dicted 30  percent  when  256-second  (64-point)  filters  were  used,  r.nd  could 
have  predicted  20  percent  with  80-second  (20-point)  iilters.  This  is,  of 
course,  only  one  example.  Each  case  shows  different  results,  with  short 
filters  sometimes  being  as  effective  as  the  longer  ones. 


5.3  HORIZONTAL  SEISMOGRAPHS 

The  discussion  of  prediction  filtering  up  to  this  point  applies  to  its  appli- 
cation to  vertical  seismographs.  As  mentioned  previously,  the  horizontal 
seismographs  at  McKinney  were  noisy  even  during  windless  periods.  A large 
number  of  single  and  multichannel  filters  computed,  and  the  results 

proved  extremely  variable  for  both  windy  and  windless  time  periods.  No 
general  conclusions  can  readily  be  deduced  from  these  experiments;  at  times 
the  single-channel  filter  was  as  effective  as  the  multichannel  filters,  while 
the  opposite  effect  was  almost  as  common.  Furthermore  there  were  numerous 
times  when  no  appreciable  improvement  was  obtained  from  any  of  the  filters. 

At  no  time  was  the  noise  level  of  the  horizontals  reduced  to  the  noise  level 
obtained  on  a horizontal  seismograph  located  at  depth.  Figures  15  and  16  show 
ordinary  and  multiple  coherences  between  the  microbarographs  and  the  north 
horizontal  seismograph  for  the  small  (100  m) , 5-element  microbarograph  array 
located  around  the  vault  at  McKinney,  Texas.  These  two  examples  are  from  windy 
periods  (approximately  6 m/sec)  when  the  coherences  obtained  were  high  and 
the  optimum  filters  were  quite  effective.  In  both  cases,  the  multiple  coherence 
and  the  highest  ordinary  coherence  of  the  5 microbarographs  is  plotted.  In  one 
case  (figure  15)  the  microbarograph  to  the  south  had  the  highest  coherence 
with  the  seismograph  while  in  the  other  case  (figure  16)  the  vault  microbaro- 
graph had  the  highest  correlation.  It  should  be  noted  that  in  both  cases  there 
was  at  least  one  microbarograph  in  this  small  array  that  showed  practically  no 
correlation  with  the  seismograph.  This  result  over  a small  array  is  not 
completely  understood,  but  it  is  usually  (though  not  always)  those  microbaro- 
graphs whose  azimuth  from  the  seismograph  is  transverse  to  the  orientation 
of  the  seismograph  that  show  low  coherences.  This  is  to  be  expected  because 
earth  tilting  in  this  direction  should  rot  affect  the  seismometer. 

Examination  of  the  figures  indicates  that  optimum  filtering  was  quite  effective 
during  these  time  periods.  However,  in  one  case  while  the  south,  north  and 
vault  microbarographs  contribute  substantially  to  the  predictability,  the 
east  and  west  microbarographs  add  little.  In  the  other  case,  the  single- 
channel optimum  filter  is  as  effective  as  the  multichannel  filters  and 
computation  of  a multichannel  filter  would  in  no  way  improve  the  results. 
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Figure  15.  Multiple  and  ordinary  coherence  for  the  north  horizontal  seismograph  and  the 
100  m,  5-element  microbarograph  array,  McKinney,  Texas. 
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6.  ATTENUATION  OF  ATMOSPHERICALLY-GENERATED  NOISE  WITH  GEOLOGY  AND  DEPTH 


During  only  moderately  turbulent  atmospheric  conditions  the  magnitude  of  the 
wind-generated  noise  can  be  as  large  as  30-60  mu  rms  in  the  vertical  direction 
and  600-1200  mu  rms  in  the  horizontal  directions. 

Experimental  data  demonstrate  that  the  effective  rigidity  of  the  rocks  at  or 
near  the  surface  of  the  earth  is  an  important  factor  in  determining  the  magni- 
tude of  the  transient  component.  As  a simple  illustration,  consider  the 
following  example  (see  section  2 for  theory  used).  Suppose  that  long-period 
seismograph  systems  are  installed  at  the  surface  of  elastic,  isotropic  half^  ’ 
spaces  whose  effective  rigidities  range  from  103  to  105  bars.  This  covers  .a 
spectrum  of  rock  types  varying  from  poorly  consolidated  sands  and  shales  to 
dense  limestones  and  ultrabasic  intrusives.  Available  experimental  data  indi- 
cates that  during  very  calm  atmospheric  conditions,  the  noise  level  in  the 
20-60  second  period  range  will  be  roughly  the  same  regardless  of  rock  type. 
However,  when  the  wind  begins  to  blow,  the  noise  levels  will  begin  to  vary. 

In  figure  17  we  display  the  anticipated  increase  in  the  20-60  second  noise 
level  produced  by  a 5 m/sec  wind  as  a function  of  increasing  rigidity.  It  will 
be  noted  that  the  noise  level  increases  in  the  horizontal  direction  are  generally 
much  larger  than  those  in  the  vertical.  The  reason  for  this  is  that  pressure 
generated  earth  tilts  cause  apparent  horizontal  earth  motion  which  is  approxi- 
mately equal  to  the  vertical  earth  motion  multipled  by  the  factor  where  g 

is  the  local  value  of  gravity,  T is  the  period  of  oscillation  and  C is  the 

convection  velocity  of  the  pressure  field  which  is  associated  with  particular 
period  (Sorrells,  1971,  Sorrells  and  Goforth,  1973).  In  our  calculations  we 
have  taken  C to  be  equal  to  the  mean  wind  speed.  Thus  in  our  case,  in  the 
20-60  second  period  range  apparent  horizontal  earth  motion  will  exceed  vertical 
earth  motion  by  a factor  which  varies  from  approximately  6 to  19  depending  upon 
the  period.  It  is  also  important  to  note  that  in  order  to  restrict  increases 
in  noise  levels  to  3 dB  or  less  in  both  the  vertical  and  horizontal  directions 
it  is  necessary  to  install  the  seismographs  in  rocks  with  effective  rigidities 
in  excess  about  3 x 105  bars.  Typical  rocks  with  rigidities  of  this  order 
would  be  dense  limestones  or  ultrabasic  intrusives  which  possess  a low  volume 
percentage  of  microfractures  and  cracks.  This  latter  requirement  is  quite 
important  since  results  obtained  by  Nur  (1971)  demonstrate  that  the  presence  of 
open  cracks  can  cause  drastic  reductions  in  the  effective  rock  rigidity. 

It  has  been  shown  by  Sorrells  (1971)  that  earth  noise  caused  by  wind-generated 
pressure  decays  fairly  rapidly  with  depth.  To  illustrate  this  effect  we  have 
calculated  the  noise  level  increases  produced  by  a 5 m/sec  wind  at  various 
depths  in  a homogeneous  and  isotropic  half  space.  The  results  are  shown  in 
figure  18.  For  computational  pruposes  we  have  assumed  an  effective  rigidity  of 
103  bars,  and  a Poissons  ratio  of  0.4.  Tiese  numbers  are  representative  of  very 
compliant  rocks.  Observe  that  in  this  case  placing  the  seismograph  at  a depth 
of  180  meters  is  approximately  equivalent  to  placing  it  at  the  surface  of  a 
half  space  characterised  by  an  effective  rigidity  of  about  3 x 104  bars. 
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Figure  17.  Relative  noise  level  during  windy  periods  vs.  rock  rigidity 
(wind  - 5 m/sec)  . 
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Figure  18.  Relative  noise  level  during  windy  periods  as  a 
function  of  depth  (wind  - 5 m/sec). 
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In  section  2,  the  theory  explaining  the  attenuation  of  atmospheric  pressure 
variations  was  briefly  discussed.  Experimental  evidence  to  confirm  the  theory 
is  still  limited,  although  the  results  show  general  agreement  between  them. 

The  main  reasons  for  the  lack  of  confirmation  are:  (1)  seismometers  have  not 
been  operated  at  a sufficient  number  of  the  shallow  depths  required  to  define 
the  amplitude-depth  curves;  and  (2)  for  horizontal  seismographs  the  environ- 
mental effects  at  the  surface  make  the  use  of  the  surface  instruments  as  a 
reference  point  extremely  difficult. 

The  main  concern  in  this  section  is  the  attenuation  of  wind-generated  noise 
with  depth.  The  effect  of  the  wind-generated  noise  decreases  rapidly  with 
depth;  in  most  experiments  the  seismographs  were  operated  at  depths  far 
beyond  those  where  the  theory  predicts  wind-generated  noise  should  be  recorded. 
The  results  of  one  of  these  experiments  (Douze  and  Sherwin,  1975)  at  Pinedale, 
Wyoming,  clearly  indicate  that  the  noise  levels  for  both  horizontals  and 
vertical  seismographs  are  approximately  the  same  for  depths  beyond  approxi- 
mately 50  m,  where  neither  wind- generated  noise  nor  surface  environmental 
effects  are  present.  Furthermore,  there  is  no  further  decrease  in  the  noise 
level  from  50  m down  to  a depth  of  914  m.  Figures  19  and  20  show  examples  of 
the  results  obtained  at  the  deep  hole  site  at  Pinedale,  Wyoming. 

Figure  19  shows  the  noise  level  of  vertical  seismographs  operating  at  depths 
of  46  m and  305  m.  The  minor  differences  at  the  longer  periods,  with  the 
shallower  instrument  somewhat  noisier,  are  not  significant  and  probably  are 
caused  by  differences  in  system  noise  levels.  Shown  in  figure  21  is  the 
coherence  between  the  vertical  seismographs;  the  rapid  decrease  in  coherence 
to  the  longer  periods  from  the  microseismic  perk  at  20-second  period  is  typical 
of  results  obtained.  The  coherert  part  must  consist  of  long-wavelength 
phenomena,  elastic  waves  in  the  earth  or  acoustic  waves  in  the  atmosphere. 
However,  at  periods  beyond  approximately  50  seconds  the  noise  consists  of 
system  noise  and  environmental  effects  that  do  not  correlate  between  seismo- 
graphs with  separations  small  compared  to  the  wavelengths  of  atmospheric 
or  seismic  noise.  Figure  20  shows  the  noise  level  of  the  horizontal  seis- 
mographs at  the  depths  of  46  and  914  m.  Note  that  the  noise  level  in  the  notch 
between  20-  and  30-second  periods  is  virtually  the  same  as  that  recorded  by  the 
vertical  seismographs  at  depths  beyond  46  m.  The  more  rapid  increase  in  noise 
level  of  the  deep-hole  seismograph  towards  the  longer  periods  is  caused  by 
convection  currents  (see  the  appendix).  At  this  depth  the  seismometer  was 
operating  in  the  water-filled  portion  of  the  casing 

At  Pinedale,  comparisons  between  the  noise  levels  at  the  surface  instrumentation 
and  the  seismographs  operating  at  shallow  depth  (46  m)  did  not  give  clear  con- 
firmation of  the  theory.  The  expected  noise  levels  were  computed  using  the 
theory  discussed  in  section  2,  the  appropriate  elastic  constants  (see 
introduction),  and  the  experimental  pressure  variations  recorded  by  the  nearby 
microbarograph. 
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Figure  19 


Power  spectra  of  the  background  noise  recorded  by  vertical 
seismographs  operating  at  depths  of  46  m (SLZ)  and  305  m 
(DLZ) . Corrected  seismograph  response;  wind  velocity 
<1.0  m/sec;  Pinedale,  Wyoming . 
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Power  spectra  of  the  noise  recorded  by  horizontal  (East) 
seismographs  at  depth  of  46  m (SLE)  and  914  m (DLE) . 
Corrected  for  seismograph  response;  wind  velocity  4.5  m/sec; 
Pinedale,  Wyoming. 
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Figure  21.  Coherence  between  vertical  seismographs  at  46 


Figures  22  and  23  show  examples  of  the  results  obtained  experimentally  together 
with  the  computed  theoretical  noise  levels  for  the  appropriate  wind  velocities. 
The  theory  for  vertical  seismographs  shown  in  figure  23  clearly  indicates  that 
a wind  of  5 m/sec  should  not  affect  the  noise  level  on  either  the  surface  or 
46  m seismograph.  Beyond  the  notch  there  is  a difference  in  the  noise  levels; 
the  power  level  of  the  surface  seismograph  increases  more  rapidly  than  that  of 
the  shallow  hole  seismograph.  However,  the  difference  is  no  larger  than  that 
often  obtained  during  windless  periods. 

For  the  horizontal  seismographs  the  data  are  more  complex.  The  surface 
horizontal  noise  level  increased  more  than  predicted  by  the  theory  in  the 
notch  (20-30  second  period) . A suspected  leak  in  the  surface  valult  may  be 
partly  responsible  for  this  behavior.  At  periods  beyond  100  seconds  the 
theoretical  atmospherically-generated  noise  level  and  the  actual  noise  level 
at  46  m are  quite  similar.  Furthermore,  the  increase  in  coherence  (see  figure 
24)  between  the  two  seismographs  at  periods  beyond  70  seconds  clearly  indicates 
that  the  two  instruments  are  responding  to  a common  input  although  the  noise 
level  at  the  surface  is  still  well  above  that  predicted  by  theory.  It  is 
possible  that  the  rigidity  at  the  surface  is  greatly  overestimated  because  of 
near  surface  fractures;  this  phenomen  could  account  for  the  discrepancy  between 
theory  and  the  experimental  evidence. 

At  the  Grand  Saline  mine  the  seismograph  at  183  m depth  recorded  no  appreciable 
amount  of  wind-generated  noise  while  the  surface  seismograph  noise  levels  in- 
creased by  large  amounts  at  this  site  in  sedimentary  material  (Sorrells  et  al., 
1971).  Thus,  no  direct  information  on  the  attenuation  of  wind-generated  noise 
was  obtained.  However,  during  the  passage  of  acoustic  vraves  from  a presumed 
atmospheric  explosion  the  results  clearly  show  that  acoustic  waves  cause  deforma 
tion  at  the  surface  and  at  183  m with  little  attenuation. 

In  conclusion,  evidence  collected  during  these  programs  cannot  conclusively 
demonstrate  the  correctness  of  the  theory  of  wind-induced  noise. 
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Figure  22.  Power  spectra  of  the  noise  recorded  by  vertical  seismographs 
at  46  m (SLZ)  and  the  surface  (LZK) . Wind  velocity  5 m/sec. 
Also  shown  is  theoretically  predicted  noise  corrected  for 
seismograph  response,  Pinedale,  Wyoming. 
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Figure  23.  Power  spectra  of  the  noise  recorded  by  horizontal  seismographs 
at  46  m and  at  the  surface.  Wind  velocity  5 m/sec.  Also 
shown  is  theoretically  predicted  noise.  Corrected  for  seismo- 
graph response,  Pinedale,  Wyoming. 
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Coherence  between  recordings  from  vertical  seismographs  at  depth  of  46 
and  at  the  surface. 


7.  CONCLUSIONS 


The  main  source  of  noise  that  affects  surface  seismographs  in  the  period 
range  of  20-100  seconds  is  the  pressure  variation  caused  by  wind-generated 
turbulence.  The  principal  proof  consists  of  the  high  coherence  between 
seismographs  and  colocated  microbarographs.  Environmental  effects  are  also 
present  and  are  especially  severe  on  long-period  horizontal  seismographs  be- 
cause of  their  tilt  sensitivity. 

The  sources  of  the  noise  during  windless  days  or  at  depth  are  not  well 
established.  Both  seismic  waves  and  atmospheric  pressure  variations  appear 
to  contribute  to  the  noise  levels.  Acoustic  waves  in  the  atmosphere  are  the 
main  contributor  and  are  the  principal  reason  for  the  seasonal  variations  in 
the  otherwise  stable  noise  spectrum. 

Of  the  two  methods  of  suppressing  wind-generated  noise,  burial  at  depth  and 
prediction  filtering,  the  former  is  the  more  effective.  At  depths  at  which 
the  wind-generated  noise  is  eliminated  (approximately  50  m in  hard  rock)  both 
horizontal  and  vertical  seismographs  have  low  noise  levels.  No  further 
reduction  is  obtained  by  installing  seismographs  at  greater  depths  than  that 
required  to  eliminate  the  short  wavelength  pressure  variations  associated  with 
the  wind.  Prediction  filters  using  colocated  microbarographs  to  predict  and 
thus  eliminate  the  wind-generated  noise  recorded  by  surface  seismographs  are 
very  effective  when  applied  to  vertical  seismographs.  For  horizontals  the 
results  are  variable  and  prediction  does  not  result  in  appreciable  noise  re- 
duction. 
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CONVECTION  IN  BOKEHOLES 


During  the  course  of  operating  long-period  seismographs  in  boreholes  it  was 
found  that  convection  currents  in  the  fluid  in  the  casing  can  seriously  affect 
the  performance  of  the  system.  The  convection  currents  are  caused  by  tempera- 
ture gradients.  Hales  (193/)  derived  a formula  for  the  instability  threshold 
for  positive  geothermal  gradients: 


(£) 


gaT 


critical 


Bvk 

gaaL 


Here 


g 

a 

T 

C 


accel.  of  gravity 
coeff.  of  expansion 
absolute  temp, 
specific  heat 


B = 216  (constant) 

v = kinematic  visosity 
k = thermal  diffusivity 
a = radius  of  hole. 


The  first  term  of  this  formula  is  the  adiabatic  gradient  of  the  fluid  (10  C/km 
for  dry  air  and  0.2°C  for  water  at  20°C)  . The  geothermal  gradient  in  the  earth 
fs  generally  greater  than  10°C/km;  hence,  stability  or  instability  will  depend 
mostly  upon  the  second  term,  which  is  strongly  dependent  on  the  hole  radius 
and  less  so  upon  the  thermal  properties  of  the  fluid.  Diment  (1967)  shows  that 
the  fluid  in  nearly  all  deep  wells  should  be  unstable.  Indeed,  he  measured 
continual  temperature  fluctuations  up  to  0.05°C  in  amplitude  and  between  1 and 
50  minutes  in  period  at  most  depths  in  such  a well  filled  with  water.  (His 
instruments  were  incapable  of  following  fluctuations  of  shorter  periods.) 
Gretener  (1967)  obtained  similar  results.  Both  of  these  workers  found  indict 
tions  that  the  convection  eddies  believed  to  exist  were  comparable  to  the  hole 
diameter  in  vertical  extent.  Both  Gretener  (1967)  and  Garland  and  Lennox  (1962) 
showed  that  the  observed  theimal  instabilities  could  be  stopped  by  reducing  the 
hole  diameter  by  means  of  a cemented-in  inner  casing  or  a loose  bundle  f t]^es- 
The  critical  diameter  corresponded  in  order  of  magnitude  to  that  predicted  by 

Hales'  formula. 

In  two  instances  the  performance  of  long-period  seismographs  was  seriously 
degraded  by  convection  currents.  The  first  case  occurred  at  the  Alaskan  ong 
Period  Array;  at  this  site  tiiaxial  long-period  seismographs  were  installed  in 

13  in.  casings. 

During  the  winter  months  the  noise  level  of  the  seismographs  increased 
appreciably.  After  much  experimentation  the  noise  was  found  to  be  caused 
by  intense  convection  currents  in  the  top  few  feet  of  the  casing.  urmg  e 
winter  a very  high  positive  temperature  gradient  can  be  expected  because  of 
the  sub-zero  weather.  A microbarograph  measuring  pressure  variations  inside 
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the  casing  recorded  pressure  variations  of  100  ubar^/Hz  at  60-second  period} 
this  is  approximately  equivalent  to  the  pressure  variations  that  occur  in  the 
atmosphere  when  the  wind  is  blowing  at  5 m/sec. 

During  operation  of  the  KS  seismographs,  one  of  the  most  prominent  features 
on  the  recordings  were  oscillations  at  periods  greater  than  60  seconds  when  the 
instrument  was  in  the  water-filled  part  of  the  hole.  Using  Hale's  (193) 
formula  the  following  critical  gradients  were  computed. 

The  effect  of  various  fluids  is  summarized  for  the  large  diameter  section  of 
the  Pinedale  well  in  the  following  table  of  calculated  critical  thermal 
gradients,  after  Hales  (1937). 


8.92  in.  diameter  borehole 


Hole  fluid:  Air  Water 

Critical  temperature 

gradient,  deg  C/meter:  0.0164  0.0003 


The  geothermal  gradient  in  the  well  is  approximately  .01°C/m,  from  a temperature 
log  table  when  the  well  was  first  drilled.  Thus,  the  geothermal  gradient  suggests 
that  the  water  column  is  unstable  while  the  casing  filled  with  air  is  stable. 
However,  the  down-hole  electronics  dissipates  heat,  thus  the  stability  of  the 
air  column  cannot  be  taken  for  granted.  The  spectra  of  the  noise  taken  for  the 
seismographs  operating  in  air-filled  casing  do  not  indicate  the  presence  of  con- 
vections of  sufficient  magnitude,  if  any,  to  effect  the  performance  of  the 
seismographs. 

Figure  zO  .In  the  body  of  the  report  shows  the  power  spectra  of  the  noise  recorded 
by  the  east  components  of  the  shallow-  and  deep-hole  seismographs.  The  deep- 
hole  instrument  was  at  914  m in  the  water-filled  portion  of  the  hole.  The 
long-period  noise  that  is  being  attributed  to  convection  currents  is  clearly 
noticeable  for  periods  greater  than  approximately  40  seconds,  and  the  power  level 
increases  rapidly  towards  the  longer  periods.  A comparison  with  the  spectra 
from  the  deep-hole  verticals  shows  that  the  long-period  oscillations  are  not 
present  on  the  vertical  component.  Therefore,  this  type  of  noise  is  almost 
certainly  caused  by  tilting  of  the  instrument  paCKage. 
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